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Abstract 

Chimera states in networks of coupled oscillators occur when some fraction 
of the oscillators synchronise with one another, while the remaining oscillators 
are incoherent. Several groups have studied chimerae in networks of identical 
oscillators, but here we study these states in heterogeneous models for which 
the natural frequencies of the oscillators are chosen from a distribution. For 
a model consisting of two subnetworks we obtain exact results by reduction 
to a finite set of differential equations, and for a network of oscillators in 
a ring we generalise known results. We find that heterogeneity can destroy 
chimerae, destroy all states except chimerae, or destabilise chimerae in Hopf 
bifurcations, depending on the form of the heterogeneity. 

Synchronisation of interacting oscillators is a problem of fundamental 
importance, with applications from Josephson junction circuits to neuro- 
science [iW, 24, 22, 27J. Since oscillators are unlikely to be identical, the 
effects of heterogeneity on their collective behaviour is of interest. One 
well-studied system of heterogeneous phase oscillators is the Kuramoto 
model [¥, 23j, for which there is global coupling. Generalisations of 
this model with nonlocal coupling [3, 21, ^ ,10, .16j, or several popula- 
tions of oscillators [l] , have shown interesting types of behaviour referred 
to as "chimera" states in which some oscillators are synchronised with 
one another while the remainder are incoherent. Even though the ef- 
fects of heterogeneity on synchronisation have been emphasised in the 
past [4*, 13, ^, 23], all chimerae have so far been studied in networks of 
identical oscillators. This raises the obvious question: do chimerae exist 
in networks of nonidentical oscillators? Here we address the question 
analytically, first using recent results to exactly derive a finite set of dif- 
ferential equations governing the dynamics of chimerae in two coupled 
networks of heterogeneous phase oscillators, and then using a similar 
idea to extend the results of Abrams and Strogatz [3j for chimera states 
in a ring of coupled oscillators. 



1 Introduction 

Networks of coupled oscillators have been studied for many years [24l [131 HI [3 
m [ini [S]. One well-known system is the Kuramoto model [1 [i [H [TOl US] of 

phase oscillators. In the last few years several authors have studied "chimera" 
states in networks of identical Kuramoto oscillators, in which some oscillators are 
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synchronised with one another while the remainder are incoherent. Much analytical 
progress has been made in the study of these states [H [31 [2T1 ^ [TOl [H] . It is very 
unlikely that any physical system being modelled by a network of coupled oscillators 
will have identical units, so the robustness of chimerae to network heterogeneity is 
naturally of interest. Certain networks of coupled oscillators are known to have 
non-generic properties |25( I26j , and it is of interest to know whether chimera states 
are generic and stable (and thus expected to be generally observed) or not. 

Here we conduct an analytical investigation into the robustness of previous- 
ly-studied chimerae with respect to heterogeneity in the intrinsic frequencies of 
oscillators. We find that chimerae are robust with respect to this type of hetero- 
geneity, and show some of the bifurcations that chimera and other states undergo as 
the oscillators in the network are made more dissimilar. Our results provide more 
evidence that the Ott-Antonsen ansatz 17j correctly describes attracting states in 
Kuramoto-type networks when the oscillators are not identical. Some of the ideas 
here have recently been used by others to study a single population of oscillators 
with a bimodal frequency distribution [14 and the periodically forced Kuramoto 
model [Ij. 

Note that the term "chimera" has been used in the past to refer to certain states 
in networks of identical oscillators [iSl [HI [2 [3] , but here we also use the term to 
describe similar states in heterogeneous networks in the obvious way. A state found 
by numerically continuing from a chimera state in a network of identical oscillators 
is also referred to as a chimera. 

In Sec. [21 we present the first model of two coupled networks, then consider its 
continuum limit and use the remarkable recent result of Ott and Antonsen [17j 
to derive three ODEs which exactly describe some of its behaviour. In Sec. [3| we 
perform a limited bifurcation analysis of these ODEs and interpret the results. In 
Sees. [4|and[5| we consider other distributions of the intrinsic oscillator frequencies, 
and generalisations, respectively. In Sec.[6|we consider oscillators on a ring. 

2 Two coupled networks 

We first consider two networks of coupled oscillators with uniform coupling between 
oscillators within each network, and a weaker coupling to those in the other network. 
Our model equations are 

■m=l j=l 

ior i — 1, . . . N and k — 1,2, where the natural frequencies uif are chosen from a 
distribution gk{ui^)- Our system is the same as that of Montbrio et al. [15]; the 
system of Abrams et al. [T] is a special case of that studied here. A similar system 
was studied by Barreto et al. [^, but their focus was the onset of synchrony, as 
was Montbrio et al's. Like Abrams et al. [1], we choose Kn = K22 = M and K12 — 
K21 — set fi + h' = I (by rescaling time if necessary), and choose fi > v. We define 
A = fi—v and (3 = -k /2 — a. Abrams et al. found that for (3 and A sufficiently small 
and positive (and all equal), both the completely synchronised state (0™ = 0^ 
for all j,i,m,k) and the chimera state (all oscillators in one population perfectly 
synchronised, all oscillators in the other population incoherent) were stable. 

We take the continuum limit of letting N ^ 00. The system is then de- 
scribed by the probability density function (PDF) /j. (a;'° , 6**^ , t) for each population 
k. We define two order parameters 

/oo /'27r 
/ exp(^0'=)/fc(c^^0^^) d^'^ (2) 
00 Jo 
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for k — 1,2. Each fk satisfies a continuity equation 

f + ^(/.^.) = (3) 

where the velocity 

- exp [t{0'' + a)]{fiZk + iyzk')}/{2i), (4) 

fc' = 3 — fc, and an overbar denotes the complex conjugate. Writing fk as a Fourier 
series in 6'^ we have 



27r 



1 + <^ ^/l„(w^^) exp {ind ~) + c.c. 



(5) 



where "c.c." denotes the complex conjugate of the previous term. Substituting (O 
into and (121), one can derive an infinite set of integro-differential equations for 
the hn [E]. However, Ott and Antonsen [17] noticed that for the special choice 

/l„(c.^^) = [a,(c.^^)]" (6) 

i.e. hn is ak raised to the nth power, all of these differential equations are actually 
the same, and we are left with a single PDE governing the dynamics of ak[ijj^ ^t): 

+ iuj^ak - (e*"/2)(^Zfc + vzu') 

+ ie-'"/2)ifiZk + i^Zk')al = (7) 

where 

Zk{t) = / ak{i^\t)gk{Lo'') dw'' (8) 



The ansatz ^ is not trivial, and the reduction from an infinite set of differential 
equations to one is remarkable. Ott and Antonsen [17 give more detail on the 
circumstances under which this ansatz is valid, and we discuss its usefulness in 
describing attracting states below. 

As is well-known (TS] [T71 [HI E] , if gk is a Lorentzian distribution the integral ^ 
can be evaluated analytically. Suppose that 

i.e. the u;'^ are from a distribution centred at ilk with half- width- at-half-maximum 
Dk- Then Zk{t) = ak{^k — "iDk^t) and evaluating ([7]) at w'^' = ilk — iDk we obtain 

'^ + {Dk+ ink)zk - (e*"/2)(^Zfe + vzy) 
at 

+ /2){^iZk + i^Zk.)zl = (10) 

i.e. a complex ODE for each k. Writing zi — rie~^'^^ and Z2 — r2e^*"^2 and defining 
(j) = — (j)2 we obtain the three real ODEs: 

dri f 1 — 

— = -Diri+i — ^-^ j [^ri cosa -f i/r2 cos (0 - a)] (11) 
dr2 / 1 — r'^ \ 

— = --02^2 + ( — ^ j [m?'2 cos a + vri cos {(f) + a)] (12) 



# ( rl + 1 



dt \ 2ri 

2r2 



[/iri sin a — vr2 sin {(j) — a)] + il2 

[/^r2 sina + i^ri sin ((/) -I- a)] (13) 
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where, without loss of generahty, we have set fli — 0. When Di = D2 = 0,2 ~ 0, 
we recover the results of Abrams et al. jj. In particular, ri = 1 {9} all equal) is 
invariant. If ri = 1, there also exists the perfect synchrony state (r2 = 1,0 = 0) 
and, depending on parameters (see Fig. 4 in [IJ) two other fixed points with r2 7^ 1 
(the chimerae). When they exist, one of these fixed points is a saddle while the other 
is either stable or unstable, changing stability via a supercritical Hopf bifurcation. 
We now proceed with a limited analysis of (fTT|l -(|13 p . 

3 Results 

3.1 Varying distribution widths, no frequency offset 

First consider the case when Di = D2 = D and fl2 = 0. The system ([TT1) -(fT3 |) 
possesses Z2 symmetry: (ri, r2, </>) — > (r2, ri, — 0). We choose A — 0.2 and (3 = 0.07 
so that for D = Q, there exist five fixed points of (fTT|) - P^ : the perfect synchrony 
state (ri, r2, 4>) = (1, 1, 0) and two chimerae (one stable and one a saddle) for which 
ri = 1 and 7-2 7^ 1 and 7^ 0, and their symmetrically related states. The fixed 
points of pT|) - p^ and their stability as a function of D are shown in Fig.[T] There 
are several interesting observations to be made here. Firstly (for these parameter 
values) , increasing the heterogeneity of the network first destabilises the symmetric 
state (ri = ^2), then restabilises it. Secondly, increasing the heterogeneity actually 
decreases the width of the angular distribution of the unsynchronised population in 
the chimera state (lower blue branch in panel (a) of Fig.[T]). 

3.2 Varying one distribution width, no frequency offset 

Now consider varying D2, while Di = D,2 ~ 0- The system no longer possesses 
any symmetry, so the effect of increasing D2 from zero on the chimerae with ri = 
l,r2 7^ 1 will be different from its effect on the chimerae with ri 7^ 0, r2 = 1. We 
choose A = 0.2 and (3 — 0.1, so that as before, when D2 = there exist five fixed 
points of (fTT|) - p3l) . Results are shown in Fig. [2l With five solutions to track, we do 
not show (f). Also, even though D2 < is not physically meaningful we plot fixed 
points for Z?2 < to show how branches of solutions are connected. 

Panel (a) in Fig. [2] shows fixed points for which ri = 1 (recall that we are 
making population 2 heterogeneous). As D2 is increased, the perfectly synchronous 
solution that exists at D2 = is destroyed in a saddle-node bifurcation, while the 
chimera with population 2 desynchronised persists. Panel (b) shows the fate of the 
two chimerae (one stable and one a saddle) for which r2 = 1 when D2 = 0. We see 
that they are both destroyed in a saddle-node bifurcation as D2 is increased. From 
this figure we see that if one population is made sufficiently heterogeneous, the only 
solution that persists is the chimera for which the oscillators in that population are 
desynchronised. Interestingly, if D2 is increased to larger values {D2 ~ 0.1), the 
state where both populations are in the "splay" state, with uniform angular density, 
i.e. ri = r2 = and (j) is no longer meaningful, becomes stable (not shown). 

3.3 Varying frequency offset Q2 

We now consider varying f22 with Di — D2 — 0. Since we are interested in states 
for which at least one of the populations is in complete synchrony we set ri = 1 and 
only consider (fT2|) - (fT3l) . The completely synchronised state, {r2,(t>) = (1,0), exists 
when n2 — 0, and as O2 is increased it persists as the fixed point {r2,(f>) — (l,(/>), 
where (p is the solution closest to zero of 2ycosasin0 = 02 ■ Numerical results 
are shown in Fig. [H We see that as O2 is increased from zero, the synchronised 
state for which r2 = 1 is destroyed in a transcritical bifurcation involving the saddle 
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Figure 1: (Colour online) Fixed points of ((TTll-lIISl) when Di = D2 = D,Vl2 = 0. 
(a): ri and r2 as a function of D. (b): (f> as a. function of D. Blue curve: stable 
chimera. Red curve: saddle chimera. Black curve: symmetric state (ri = r2). 
Solid lines indicate stable solutions, dashed lines unstable. Other parameters: A = 
0.2, /3 = 0.07. 
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Figure 2: (Colour online) Fixed points of ([Il])-(IT3l) when Di = ^2 = 0. (a): r2 
as a function of D2 when ri = 1. (b): ri (red) and r2 (blue) as a function of 
D2- Solid lines indicate stable solutions, dashed lines unstable. Other parameters: 
A = 0.2, /3 = 0.1. 
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Figure 3: (Colour online) r2 as a function of when ri = 1. Solid lines indicate 
stable fixed points of (fTT1) - P^ . dashed lines unstable. Circles are the maximum 
and minimum of r2 during stable periodic oscillations. Note that the branch with 
r2 > 1 is not physically meaningful. Other parameters: A = 0.2,(3 = 0.1. 



chimera, while the stable chimera undergoes a supercritical Hopf bifurcation, leading 
to oscillations in r2 and 4>. However, decreasing Q2 from zero causes destruction 
of the stable chimera in a saddle-node bifurcation with the saddle chimera. For 
these parameter values, one can see that the stable chimera is much more robust to 
speeding up the asynchronous oscillators, as opposed to a slowing them down. 

3.4 Discussion 

Our bifurcation analysis in this section has found all four codimension-1 bifurcations 
of ODEs (saddle-node, pitchfork, transcritical and Hopf) and a two-parameter study 
is likely to find higher codimension bifurcations. 

Pikovsky and Rosenblum [18] recently studied ([T|) with identical (i.e. the 
system of Abrams et al. [T] , and our system when Di = D2 = O2 = 0) as a special 
case and found that the ansatz ([6]) did not completely describe the possible dynamics 
of this system. However, Pikovsky and Rosenblum [18] and other authors [M] found 
that when the oscillators are non-identical, this ansatz does successfully allow one 
to describe attracting states. We also find this behaviour here: extensive numerical 
simulations show that the stable states shown in Figs. [1] and O for D,D2 > are 
attracting, and that the angular distributions of these stable states are given by (O- 
([6]), even if the initial distributions are not. However, the same cannot be said for the 
results in Fig. |3l for which oscillators within each of the two networks are identical. 
This figure correctly predicts the dynamics if the initial angular distribution is given 
by dSJ-jH]) but other initial conditions give solutions not described by Fig. [3] (not 
shown) . This relationship between initial conditions and dynamics when oscillators 
within each network are identical was also noticed by Montbrio et al. [15]. The 
results in Fig. [3] are likely to be a subset of those that could be found using the 
approach in Ref . [18] . 

In related results, Montbrio et al. [15] fixed ^2^0 and varied both a and ly 
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Figure 4: (Colour online) ri and r2 fitted to simulations of H]), where all ojf are 
chosen from a normal distribution of mean zero and standard deviation a. Blue 
circles joined by a line: stable chimera. Black crosses joined by a line: stable 
symmetric state (ri = r2). Red dashed line: presumed unstable symmetric state. 
Compare with Fig. [TJ See the text for details on the fitting. Other parameters: 
^ = 0.2, /3 = 0.07, iV= 1000. 



and found chimera states, both for homogeneous and heterogeneous networks. 

4 Other distributions 

Now we consider the effects of choosing the uj^ from distributions other than the 
Lorentzian, first numerically and then analytically. 



4.1 Gaussian distribution: numerical simulations 

Figure [3] shows the results of fitting the time-dependent PDF 



fkiO,t) = ^ 

2,71 



l-ri 



27r[l - 2rfe cos (4>k 



(14) 



to each population in simulations of ([T]) after transients, where all u)^ are chosen 
from a normal distribution of mean zero and standard deviation a. We found that 
both ri and r2 tended to constant values, as did 02 — 4>i (not shown). Only stable 
states are shown in Fig.[H but the results are compatible with those shown in Fig.[Tl 
suggesting that there is nothing special about the Lorentzian distribution, as has 
been noted by others [14j [7] . The unstable states could presumably be found using 
the "equation-free" method [13l [H [12] of analysing low-dimensional descriptions of 
high-dimensional systems, under the assumption that these states are also exactly 
described by the variables r, (p for each population. 



8 



4.2 Another distribution 

As an alternative [IJ, we suppose that the lo^ are chosen from the distribution 



^/2s 



1 



(15) 



which has mean zero (for simpHcity) and variance s\. This gk{^) has poles at 
u) — Sfc(±l ± *)/V2, and the integral ([S]) gives 



Zk 



l + i 



1 - i 



■,+ 



(16) 



where satisfy 



dt 



Sk 



(1 ± i)/V2] z± - {e-'" /2){^xzk + z^ZfeO 
+ (e'"/2)(/^Zfe + i.ZfcO(^fc)' = 



(17) 



Thus we have four coupled complex ODEs, rather than the two ((TO)) . This can 
clearly be generalised to other distributions which are rational functions of uj. 



5 Generalisations 



We now briefly mention several generalisations of the results above. Suppose that 
the system ([T]) is periodically forced, i.e. we add the term sin {nt-e^) to (P), 
as done recently for a single population [7 . Going to a coordinate frame rotating 
with angular frequency 51, we find that (jlOp is replaced by 



dzk 
dt 



{Dk+iiVLk 



Afc + e'^iiiZk + vzk') 



n))zk 

Ak + e~"^{fiZk + i^Zk') 



zl^O 



The system is no longer invariant under a translation of time, so during the deriva- 
tion of ODEs like PT|) - P^ we find that we need both 0i and 7 not just their 
difference. 

Another possibility is that there is a uniform delay of r between the two pop- 
ulations, but zero delay within them, i.e. we replace 6*"' in ^ by 6'™(t — t) when 
m = k' . The effect of this is to replace Zk'{t) in (fTO|) by Zk'{t — T) and the equivalent 
of (|li p - (fT3|l is now the four delay differential equations: 



dri 
~dt 



dr2 
li ^ 

d<p2 
dt 



#1 
dt 

D2r2 - 



2ri 



2r2 



[/iri cos a 


\vr2it- 


t) cos — 




r) 


-a}] 


= 


[/iri sin a 


- vr2{t - 


r) sin {(/>! — 


Mt- 


r) 


-a}] 


= 


[/i7'2 cos a 


+ vrx{t ~ 


r) cos {4>2 — 


Mt~ 


r) 


-a}] 


= 


\[ir2 sin a 


- vri {t - 


t) sin {02 — 


Mt- 


r) 


-a}] 


= 



If all of the variables on the RHS of ([T]) were delayed by r, we could define = (t>i—<p2 
as before, and derive three coupled DDEs rather than four above. The analysis of 
the equations in this section remains an open problem. 
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6 Oscillators on a ring 



We now consider a ring of oscillators, with non-local coupling between them. The 
original presentation of chimerae was in such a system, with identical oscillators [31 
[21 Hn] • The chimera state for this system consists of oscillators on one part of the 
ring being synchronised, while over the remainder of the ring they are incoherent. 
Simulations reported in [3] indicate that these states are also robust with respect to 
perturbations of the oscillators' natural frequencies, and we now show how to use 
the ideas above to investigate this analytically. 

Consider the model consisting of oscillators on a ring studied in refs. [ini[5]i but 
include heterogeneity in the intrinsic frequencies of the oscillators. The system is 

-^=^^-^Z^G'( jcos{(P^-(j)J-p) (18) 

for i — 1, . . . iV, where the natural frequencies uji are chosen from a distribution 
g^Lij). The coupling function G is periodic with period 27r. Equation ()18|) is the 
discrete version of 

^=uj-J^ G{x-y)cos[<j>{x,t)-<j>{y,t)-(3]dy (19) 

which for constant uj is the same as that studied by [10[ [3]. The analysis for a 
heterogeneous network is very similar to that for a network of identical oscillators, 
so we skip many of the details here and refer the reader to [3] . The main difference is 
that UJ is now a variable, and certain quantities now have to be replaced by integrals 
over Lu, weighted by g{Lu). 

6.1 Analysis 

First we go to a rotating reference frame with angular speed fi, i.e. let 9 = 4> — fit. 
Then we define an order parameter 

R{x, t)e^^("'*) = / G{x- 2/)e'^(2/,t) 
Jo 

so that can be written 

f)f) 

— = uj-n-Rcoa{e-e~(3) (20) 

We now look for stationary states, so that R and O are independent of t. At position 
y, if R{y) > |a; — flj, then the oscillators will move to the stable fixed point 9* , which 
is given by the solution of 

u;-n = Rcoa [9* - - P] 

For those drifting oscillators at y with i?(y) < jw — n|, we replace e^^^"^^ in the order 
parameter definition with its average over 9 jSlllOj. but now weighted by g{uj) (over 
the appropriate range of uj). After some calculation the result is that at stationarity 
we have 

i?(a;)e*®(^) = 
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In some cases this double integral can be exactly evaluated. We follow [7 and 
suppose that G(x) = (1 + ylcosx)/(27r), so that G{x — y) — {1 + Acosxcosy + 
Asinxsin2/)/(27r). Let us define 



Hy) 



g(Lj)duj 



Thus under the assumption that R and O are even (which can be shown to be 
self-consistent) 



R{x)e 



i0{x) 



a cos X 



where 



and 



^^|3 /.27r 



2tt 



h{y)dy 



27r 



•'^'^y^h{y) cosy dy 



(23) 



(24) 



Since pT|) is unchanged by the shift Q{x) Q{x) + Oq, we can take c to be real. 
To write the right hand sides of ([23]) and (|24p in terms of a and c, note that 



i?2(y) = \R(y)e''^iy)] [i?(y)e-*®(^) 



2cRe(a) cos y + |ap cos^ y 



and 



^^^y^h{y) = R{y)e^iy)!^ 

R{y) 



(c + a cosy) 
1 



g{uj)duj 



3 V ^'(y^ 



c + a cosy 

where overbar indicates complex conjugate. Thus we have 



„j/3 r27r 



1 



2tt Jq c + a cos y 



and 



where 



2tt 



cosy 



a cosy J _ 



f{^iy)9{<^)duJ dy 

1 

O 

I{^,y)9{^)d^ dy 



(25) 
(26) 



/(w, y) = uj — ^ — — f7)2 — c? — 2cRe(a) cos y — |ap cos^ y 



Taking the real and imaginary parts of (j25[) and (|26p we obtain four real equations 
for the four real unknowns c, Re(a), Im(a) and f2. 
As in Sec. [21 suppose that 



s(^) 



D2 



1 

27ri 



1 



1 



— iD uj + iZ? 



(27) 



i.e. the u>i are from a distribution centred at zero with half- width- at-half-maximum 
D. (There is no loss of generality by assuming that the distribution is centred at 
zero, since if this wasn't so, the effect would just be to add a constant to Q..) Then 
for any function F{ijj) analytic in the lower half of the complex to plane. 



F{uj)g{uj)duj = F[-iD) 
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p -a 

Figure 5: Solutions of (HH) and ^ when D = 0. Left: Re(a) versus ^. Right: 
Re(a) versus — fi, for A = 0.95. 



and thus and (pS)) become 

_ -e*'^ /■^'^ 17 + zD + ^(17 + iD)^ - c2 - 2cRe(a) cos y - |a|2 cos^ y 



2tt Jq c + a cos y 

and 



(28) 



/-Sir (0 + 21)+ ^(il + - c2 - 2cRe(a) cosy - |aP cos^ y \ cosy 

dy 



27r 7o c + a cos y 

(29) 

Note that by setting 13 = in ([28|) and (f29|) we recover the results of Abrams and 
Strogatz [3]. 



6.2 Results 

We now show results of following solutions of (pS)) and using 13 and /3 as 

bifurcation parameters. It is known for identical oscillators (i.e. D — Q) that for 
fixed A > Q chimerae exist for a range < /3 < /?*, and that (3* is an increasing 
function of A (see Fig. 8 in [3]). Here we set A — 0.95. Results for Z3 = are shown 
in Fig. m four types of solution are shown. 

Blue crosses indicate the modulated drift state which occurs for /3 = 0, Im(a) — 
0. In this state none of the oscillators have synchronised with one another. The 
green solid line indicates the stable chimera, for which some of the oscillators are 
synchronised with one another while the remainder drift. The red dashed curve is 
the unstable chimera. The black dash-dot line represents the uniform drift state 
for which (3 = a = 0. Along this line (which has collapsed to a point in the left 
panel of Fig. [5]) all oscillators are synchronised, 9i = 9j,yi,j. If A was decreased, 
the saddle-node bifurcation seen in the left panel of Fig. [5] would move to a lower 
value of f3. (Note that stability of solutions is inferred, as it was by Abrams and 
Strogatz ;3^. All we have in (I28p and (f29|) are algebraic equations governing the 
steady states, with no dynamics.) 

A similar picture when D = 0.01 is shown in Fig. [51 with the same conventions. 
We see that the modulated drift state (which no longer has Im(a) = 0) has moved 
away from /3 = 0, as has the uniform drift state (a = 0). Although we have 
not indicated it in Fig. O for D = and < /3, the synchronised state (with 
a = 0) is stable. However for < Z3 there is now a range of P (approximately 
0.03 < f3 < 0.165 when D — 0.01) for which the synchronised state is unstable. 
(This is the extent of the dash-dotted line in the left panel of Fig. [51) For /3 outside 
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Figure 6: Solutions of ((28]) and ^ when D = 0.01. Left: Re(a) versus /3. Right: 
Re(a) versus — fi, for A = 0.95. 



this range, i.e. for /3 small enough or large enough, the synchronised state remains 
stable. 

Comparing the left panels of Fig. [5] and Fig. [6] we see that the case of identical 
oscillators {D — 0) is degenerate in the sense that both the modulated and uniform 
drift states arc "hidden" at /? = 0, but both occur over finite intervals of /3 when 
0<D. 

Fig. [7] shows the results of following the pitchfork and saddle-node bifurcations 
seen in Fig. [6] as D is varied. Note that the two pitchfork bifurcations emanate 
from {D,f3) = (0,0). The rightmost pitchfork bifurcation changes from sub- to 
super-critical at the termination of the curve of saddle-node bifurcations. For ap- 
proximately 0.047 < D < 0.058 there is no bistability; instead there is a range of 
f3 values for which only the chimera state is stable. Outside this range only the 
synchronous state is stable. For D greater than about 0.058 there do not exist any 
chimera states. 

From Fig. [7] we see that with j3 small and fixed, increasing D first destabilises 
and then restabilises the synchronous state. This is the same behaviour as observed 
in Fig. [l]for the model H]), and is demonstrated in Fig. [5] where we fix /? = 0.15 
and successively increase D from to 0.02 to 0.06. For D = the synchronised 
state is stable (bottom left panel). At t ~ 100 a chimera forms, with the center of 
the unsynchronised cluster at i « 200 and the center of the synchronised cluster at 
i « 700 (bottom middle panel). Once D is increased above the upper blue curve in 
Fig. [7] only the (noisy) synchronised state is stable (bottom right panel). 



6.3 Generalisations 

As shown in Sec. 14. 2[ if 



V2D^ 



we could repeat the analysis of (|19p . obtaining equations similar to, though more 
complicated than, (pS)) and still with the unknowns c, a and Vl. Indeed, for any 
distribution g the double integrals (|25p and could be evaluated numerically. 

The form of (|22p is a direct result of choosing G to have one Fourier mode. 
If more modes were used in G (if, for example, we were approximating a given 
coupling function with a finite Fourier series), eqn. ()22[1 would have more terms and 
thus more coefficients to be found. 
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Figure 7: Curves of pitchfork (solid) and saddle-node (dashed) bifurcations of solu- 
tions of (|^ and Fig. [5] corresponds to a horizontal "slice" through this figure 
at D = 0.01. A = 0.95. 




500 1000 500 1000 500 1000 
Index Index Index 



Figure 8: Results of a simulation of (fT5|) for which D is switched from to 0.02 at 
t = 50, and then increased to 0.06 at t = 300. Top: sini^^. Bottom row: average 
of d(j)i/dt as a function of index i over the time intervals [0,50] (left), [150,300] 
(middle) and [300,450] (right). TV = 1000,/? = 0.15, A = 0.95. 
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6.4 Chimera states and "bumps 



Chimera states as studied in this section are very similar to "bump" states which 
have been studied in computational neuroscience modelling [121 lllj . The main 
difference between bump states in neural models and the chimera states studied 
here is that in neural models, the uncoupled unit is a model neuron which — as an 
input current is increased — starts to fire periodically once the current has passed 
a threshold |11| . whereas in the chimera states studied here the uncoupled unit is 
a phase oscillator with uniform angular velocity. In neural models, a bump is a 
self-consistent state for which some neurons receive subthreshold input (and are 
thus quiescent) while other receive superthreshold input (and are thus firing). It 
is the coupling of phase oscillators through a sinusoidal function of the phase itself 
(as in (j20p ) which allows some oscillators to lock and rotate at a uniform frequency 
(which can be set to zero by moving to a rotating coordinate frame), while others 
drift, and thus a chimera can form. 

To further emphasise the similarity, compare Fig. 4 in [11] with the insets in 
Fig. 12 in [3], and with the middle plot in the bottom row of Fig. [8] (keeping in 
mind that this is a disordered system). The bifurcations of bumps are also very 
similar to those of chimerae on a ring — they both typically appear as unstable 
states bifurcating from a spatially-uniform state; compare Fig. 10 in with Fig. 12 
ini. 

7 Summary 

We have considered the effects of heterogeneity in the intrinsic frequencies of os- 
cillators on chimera states in Kuramoto-type networks of coupled phase oscilla- 
tors. Previous authors had only considered these states in networks of identical 
oscillators [TJ [21 [31 UHl HEl HI] • By assuming a Lorentzian distribution of intrinsic 
frequencies we have generalised the results of Abrams et al. Ij and Abrams and 
Strogatz [21[3], obtaining similar equations to them, but with an extra parameter, 
viz. the width of the Lorentzian distribution. All of our results show that chimerae 
are robust — within limits — to heterogeneity in their intrinsic frequencies, and 
we have shown some of the interesting bifurcations that can be induced by such 
heterogeneity. 

Importantly, in light of the recent results of Pikovsky and Rosenblum [TB] re- 
garding the validity of the Ott-Antonsen ansatz ^ used in this paper, our numerical 
results in Sec. [3] support the observation by Martens et al. [Ml that this ansatz can 
be used to study all attractors of a Kuramoto-type system whenever the oscillators 
have randomly distributed frequencies. 

The results presented here rely on the form of the equations studied. In par- 
ticular, the results in Sees. [210 rely on the remarkable recent results of Ott and 
Antonsen 171 showing that the infinite network can be exactly described by a finite 
number of ODEs, although not necessarily completely [18]. Similarly, the analysis 
in Sec. [6] depended on the form of the coupling in (fTS]) . through a trigonometric 
function of phase differences. The challenge remains to discover similar results for 
oscillators not described by a single variable, and not coupled in this way. 

Acknowledgements: I thank the referees for their very helpful comments. 
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